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Abstract 



We show numericahy that the nature of the doping induced metal-insulator 
transition in the two-dimensional Hubbard model is radically altered by the 
inclusion of a term, W, which depends upon a square of a single-particle 
nearest-neighbor hopping. This result is reached by computing the localiza- 
tion length, in the insulating state. At finite values of W we find results 
consistent with S,i ^ — /ic|~^^^ where //c is the critical chemical potential. 
In contrast, ~ 1/^ ~ /^c| "'^''^ for the Hubbard model. At finite values of W, 
the presented numerical results imply that doping the antiferromagnetic Mott 
insulator leads to a d^2_y2 superconductor. 
PACS numbers: 71.27.-ha, 71.30.-Fh, 71.10.-hx 

I. INTRODUCTION 

The doping induced metal-insulator transition in the two-dimensional Hubbard model 
is anomalous in the sense that it cannot be understood in terms of a generic band metal- 
insulator transition |l|-0. This statement follows from numerical work at zero- and finite- 
temperatures on lattice sizes up to 16 x 16. Zero-temperature quantum Monte-Carlo (QMC) 
calculations of the charge susceptibility, Xc = dn/dfi, have shown that it behaves as Xc ~ 
l/U — yUc|~^/^ in the vicinity of the critical chemical potential /ic at which the transition occurs 
IQ]. This result is consistent with t — J model calculations In principle, the square 

root singularity in Xc can be reproduced within a mean-field spin density wave (SDW) 
approximation provided that perfect nesting is present and that antiferromagnetic order is 
retained upon shifting the chemical potential in the rigid band. When perfect nesting is 
absent the SDW approximation yields, as in the case of the generic band metal insulator 
transition, a chemical potential independent charge susceptibility. This stands in contrast 
to QMC results, which for a small violation of perfect nesting still show a square root 
singularity in Xc as the critical chemical potential is approached 0]. The second piece of 
information signaling the anomalous nature of the metal insulator transition, is obtained 
via the numerical calculation of the localization length, in the insulating state. For the 
two-dimensional Hubbard model one obtains ~ |/i — /id"^''^ 0- In contrast the SDW 
approximation as well as the generic band metal-insulator transition yield a square root 
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singularity for ^i. Finally, a QMC calculation of the high frequency Hall coefficient, shows 
no divergence of this quantity in the vicinity of the metal- insulator transition |^ . As opposed 
to the case of the generic band metal insulator transition, it's magnitude remains small. 

In a single-particle picture, those numerical observations may be understood by the oc- 
currence of flat bands around the (ivr, 0), (0, ±7r) points in the Brillouin zone. More precisely, 
at the critical chemical potential, the exponent obtained numerically from the localization 
length, imposes a dispersion relation around those points. This flat dispersion rela- 
tion required for the understanding of numerical data, has important implications. It must 
originate from a large wavenumber dependent renormalization of the charge excitations due 
to correlation effects since the ordinary van-Hove singularity does not generate such a pro- 
nounced flatness. Because of this flatness around the (ivr, 0) and (0, ivr) points, doped 
holes predominantly occupy those regions. Thus, at low doping the the regions around the 
(±7r, 0) and (0, ivr) points in the Brillouin zone govern the low energy physics. Such strong 
wavenumber dependent renormalization may in part be caused by strong scattering of the 
quasiparticle between (vr, 0) and (0, vr) by the antiferromagnetic fluctuation carrying mo- 
mentum transfer (tt, tt) as well as by Umklapp scattering [|T^. The flat dispersion thus 
generated further enhances this scattering channel in a consistent fashion. The experimental 
observation of flat bands, in high-Tc cuprates is documented in |TT|-|T3|. 

The occurrence of flat bands suppresses the kinetic energy in the metallic state close to 
the metal-insulator transition. This places the metallic state as well as the nature of the 
metal-insulator transition close to potential instabilities. In this article, we show that the 
inclusion of a term W which depends upon the square of a nearest neighbor hopping changes 
radically the above described properties of the metal-insulator transition. We argue that 
the low energy physics of the added term, W, is partially contained in the three site terms 
T^ obtained in a strong coupling expansion (second order perturbation in the hopping) of 



the Hubbard model. 

Our starting point is the Hubbard model on a square lattice: 

Htu = -^E^?+ f^EKr - - ^) (1) 

i i 

with the hopping kinetic energy 

Kj=y f ct c^, ? +(i t ci)- (2) 



Here, cl^ (Ci-o-) creates (annihilates) an electron with z-component of spin a on site 

ci^c^^, and 6 = ±ax, ±ay where a^, ay are the lattice constants. The energy will be measured 
in units of t. The interaction we add, is given by: 

Hw = -WY.Kl (3) 

i 

One motivation of considering this form is that at half-band filling (/i = 0) the Hamiltonian: 

Htuw = Htu + Hw (4) 
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is particle-hole invariant, and the IV— term has a simple Hubbard Stratonovitch transfor- 
mation. This has for important consequence that the sign problem may be avoided in QMC 
approaches at half-filling. The added 14^— term may have various origins, as pointed out 
in Ref. fl^. One possible route is to consider the nature of the metallic state near the 
Mott transition described above. The Mott insulating state itself is characterized by the 
irrelevance of single-particle processes due to the opening of a quasiparticle gap. In the 
insulating state, the occurrence of long-range antiferromagnetic order shows the relevance 
of two-particle processes in the particle-hole channel. In two-dimensions, the flatness of 
the charge excitations at small doping suppress the relevance of single-particle processes 
and enhances the importance of two-particle processes. In contrast to the Mott insulating 
state, particle-particle processes now contribute to the two-particle channel. The form 
contains such generic terms of two-particle processes. 

Here, we show the relationship between the W-term and the three site terms obtained 
in a strong coupling expansion of the Hubbard model. The ly-term may be written as 

Hw = Hy) + + + (5) 

with 

i/W = -SWN (6) 
where N denotes the number of sites. 

and 



Hli^ = -WT.Hs'^ls + ^ls'^,:}- (9) 

Here, , = clA . , tJ = ct cl ^, , tJ „ = (clA ^) /V2, and a1 = 

) /V2. In the presence of the Hubbard U with U > W, only H^y should 

(2) 

dominate low energy physics. H^r may be neglected since it involves processes containing 

(3) 

double occupied sites costing an energy U. Hy^ describes triplet pair-hopping. The fact that 
this term comes with a positive sign, and that local triplet configurations are not favored by 
the superexchange generated by the Hubbard [/, gives us arguments to conclude that it will 
not affect the low energy physics. Thus, is the only relevant term. It may be rewritten 
as: 



i/g) =2WY. {S. ■ S.^g - 1 (n. - 1) {n,^^ 



;,<5 

^ V n - t t 
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The first term in the above equation is obtained by setting 5 = S' in and corresponds 
to a superexchange. Here S'j denotes the spin operator on site i. The second term of the 
above equation corresponds precisely to the three site term |15|. The third term is just the 
particle hole transform of the three site term. 

The derivation of the relevant terms in via a second order expansion in t/U is along 
the same line of idea as the above mentioned necessity to include two-particle processes. 
However, the strong coupling expansion by itself contains neither the physics of fiat-band 
formation nor the unusual properties of the metal-insulator transition. In that sense, the 
VT-term emerges as relevant processes only after complicated many-body processes of the 
flat-band formation and the second order result provides only a starting point. 

In the half-filled case, where QMC simulations are free of the sign-problem, the model 
(I) has been studied as a function of W/t and at fixed value of the Hubbard repulsion, 
U/t = 4 [14|. It has been found that this model shows a quantum transition at Wc/t ~ 0.3 
between an antiferromagnetic Mott insulator and d^^-y^ superconductor. Here, we fix the 
values of U and W to U/t = i and W/t < Wc/t and study the model as a function of 
chemical potential. The organization of the article is the following. In the next section, 
we study the metal-insulator transition by approaching the critical point from the insulator 
side. We review the definition of the localization length and compute it for the parameter 
set U/t = 4 and W/t = 0.05 and show that it diverges as |/i — nd'^^"^- This is the central 
result of the section and stands in contrast to the W/t = results: ~ |/i — jJc\~^^^- We 
discuss the relation between the transition in correlation length exponent between W/t = 
and W/t > 0.05 in terms of the single-particle density of states. In Sec. ( PD we consider 
finite doping. Here, we are confronted to a sign problem, which makes it exponentially hard 
to reach low temperatures and large lattice sizes. We show that at fixed hole density and 
Hubbard repulsion, the c?— wave pair-field correlations increase substantially with growing 
values of W/t. In contrast no such increase is observed in the extended s— wave channel. 
We equally study the evolution of n{fi) as a function of W/t. The last section is devoted to 
discussions and conclusions. 

From the technical point of view, our T = data is obtained with the projector QMC 



(PQMC) algorithm [16,171 . A numerically stable calculation of imaginary time dependent 



Green functions within this algorithm is described in Ref. |]TB[. Finite temperature grand 



canonical QMC algorithms were equally used [19,20|. The generalization of those algo- 



rithms to efficiently incorporate the H^— term is reviewed in |TJ]. Dynamical quantities were 
obtained with the use of the Maximum Entropy method reviewed in pl|,|22[]. 



II. INSULATOR-METAL TRANSITION AS A FUNCTION OF CHEMICAL 

POTENTIAL 

We start by studying the metal-insulator transition by approaching the critical point 
from the insulating side. This method has the great advantage that it allows us to avoid 
the sign problem in the QMC method. This method has already been applied to the case 
of the single band two-dimensional Hubbard model 0. There are implicit assumptions in 
this approach, i) The insulating state, is defined by the vanishing of the Drude weight or 
charge stiffness, D. The models we consider are particle- hole symmetric at yU = and thus, 
the critical chemical potential, fic, at which the transition occurs is defined by: 
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(10) 



For the range of chemical potentials for which the Drude weight vanishes, we will assume 
that the charge susceptibihty equally vanishes: 

X.(/.) = ^ = for < (11) 

This excludes the possibility of localized states within the charge gap as defined by the 
vanishing of the Drude weight. The occurrence of localized states within the charge gap 
could happen in the presence of disorder, ii) Our approach is based on the knowledge of 
the single-particle Green function. Thus, we have to assume that the quasiparticle gap, 
Agp = lim^v^oo Eq^^ — Eq , satisfies: 

Agp = /i, (12) 

Here, denotes the number of sites of the considered square lattice and Eq is the ground 
state energy in the A^-particle Hilbert space. In the case of an insulator to s-wave supercon- 
ductor transition, the difference between Agp and fic is of the order of the superconducting 
gap at the critical chemical potential. In this case, the insulator-superconductor transition, 
should be studied with the two-particle Green function. In the case of a transition between 
an insulator and rf— wave superconductor the nodes in the superconducting gap lead to 
Agp = He- As we will see, in the model of Eq. the d— wave pair-field correlations domi- 
nate over the s-wave pair-field correlations at finite doping so that our approach is justified. 
In the case of a insulator-metal transition where the metallic state is described by a Fermi 
liquid, the relation fic = Agp clearly holds. 

With the above assumptions, the localization length we consider is computed as follows. 
The zero-temperature Green function, defined by 

G{f,uj) =tJ2 f dte^-*(^o|rc,,r(t)c^g(0)|^o), with c^,r<t) = e^^*c,,,-€-*^*, (13) 

is real for values of \uj\ < fic and satisfies: 

G{r,uj = ^) exp{-\r\/^i) where 6 ~ |yU - yUcI"'' (14) 

Here, |\l/o) denotes the ground state of the Hamiltonian H = Htuw (see Eq. (|D) at half-band 
filling. The correlation length exponent v gives us valuable information on the nature of the 
insulator-metal transition, may be interpreted as the localization length of an impurity 
state at energy |a;| < Hc- 

In order to obtain the correlation length exponent, we need an accurate evaluation of 
the critical chemical potential, as well as the the single-particle Green function G{r, uo) for 
\uj\ < fi. Both those results may be obtained by using the PQMC algorithm described in 
to compute the imaginary time single-particle Green function 

ri~ , ni ^'^<>\^'<''''(MJ'^«) „, , (1'olE,cL,,,(-T)ca_,|1-o> 
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where Cf,a{T) = e'^^Cf^a^''^^ . Here, we use the same notation as in Eq. ([I3D and 6(t) 
denotes the Heaviside function. To obtain an estimate of the quasiparticle gap, Agp, we fit 
the tail of G{f = 0, r) to a single exponential form: exp(— rAgp). Alternatively, one may 

compute Agp by fitting the tail of 'Ek,e(k)=o ^(^' ^)/ '^k,e(k)=o exp(-rAgp). Here, 

e{k) = —2t{cos{kx) + cos{ky)) and e{k) = determines the non-interacting Fermi line in 
half-filled case. The so obtained values of are presented in Fig. |1] for W/t = 0,0.05 and 
0.15. As argued above, /ic = Agp in the model of Eq. (|). From Fig. |l]one notices that size 
effects become increasingly important as W/t grows. 

For |a;| < fic, we can efficiently calculate the single-particle zero-temperature Green 
function G{f, u) through the relation 

/oo 
drG(f, r)e™. (16) 
-oo 

The Green function (^(r*, r) is computed at half-band filling where the sign problem is not 
present and the statistical uncertainty does not grow exponentially with lattice size. How- 
ever, since we are multiplying the QMC data by the factor e™, the statistical uncertainty 
will grow exponentially with increasing values of r for tu > 0. For each lattice size L, we 
have considered the largest distance R = (L/2,L/2). For this distance, G{R,t) is plotted 
in Fig. ^ Due to particle-hole symmetry at /z = 0, G{R, r) = —G{R, — r). 

We carry out the imaginary time integration involved in Eq. (|l^) in the following way. 
For ;ar(?e positive values of r, G{R, r) ~ exp(-rA,p) E.(*^|Ttc(j,.T^|*r')«+'l4,.l'I'o^) 
Here, T'Ic^JTj^ = Cj^^. The Hamiltonian commutes with T^. Since we are working with 
periodic boundary conditions, T| = 1. Thus, T^I^E'o >= ±|^o >> and for large values of r, 
G{R,t) ~ riG{0,T) with rj = ±1. We can fit the tail of G(0, r) to the form aexp(— rA^p) 
and use this form to estimate the tail of G{R,t). The advantage lies in the fact that it is 
easier and more precise to extract the tail from G{0, r) than from G{R, r). We now estimate 
G{R, uj) with 

drG(^, r)e™ + r/a / dre^^'^-^"^^ - r]a dre"('^+^'''\ (17) 

-Tl Jti J —oo 

Clearly, ri has to be chosen large enough such that the left hand side of the above equation is 
Tl independent. Numerically, we choose ri such that for values for r > ri, G(-R, r) = r]G{0, r) 
within our numerical accuracy. 

The so obtained values of G{R, u = jj) are plotted in Fig. ^ at U /t = A and W/t = 0.05. 
for L = 6 to L = 12 lattices. One may fit the data to the form exp(— |i?|/^;) to obtain the 
localization length The so obtained locahzation length is plotted versus |/i — /id on a 
log-log scale in Fig. ^. For comparison, we have included our estimate of the localization 
length at U/t = 4 and W/t = 0. The central result of this section is the fact that the 
numerical data are consistent with 

~ l/i - jjcl'^/^ at U/t = 4, W/t = 
^/ ~ l/i - yUcr^/^ at u/t = 4, W/t = 0.05 (18) 

We also note that the localization length grows at fixed fi with increasing W/t. At W/t = 
0.15 growing size effects rendered the determination of the localization length difficult for 
the studied system sizes. 
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As mentioned in the introduction, the correlation length exponent u = 1/4 at W/t = 
may be understood in a single-particle picture by the existence of a \k\^ dispersion relation 
around the (ivr, 0) (0, ivr) points in the Brillouin zone. This very flat dispersion relation 
should lead to a high density of states, A^(a;) = J2k^T^G{k,u), for frequencies close to the 
quasiparticle gap. In comparison, at finite values of W/t where we obtain the correlation 
length exponent u = 1/2, we expect no anomalous enhancement of A^(ti;) at frequencies 
equal to the quasiparticle gap. In fact, for the generic metal-insulator transition for two- 
dimensional systems [^, for which one obtains the exponent u = 1/2, A^(ci;) remains a finite 
constant until it jumps to zero at the band-edge. In order to check this idea, we have 
computed the one-electron density of states: A^(a;) at half-band filling. At zero-temperature 
and Tt > 0, N{uj) is related to the imaginary time QMC data via 

G{r = 0,r) = - /"du;e"™Ar(u;). (19) 

71 Jo 

N{uj) is extracted from the above equation with the use of the Maximum-Entropy method. 
Fig. 1^ plots N{uj) for lattice sizes ranging from L = 6toL = 16atf//t = 4 and for both 
W/t = ( Figs. |f to |j ) and W/t = 0.15 ( Figs. |a to |8 ). Figs. |k to |n plots the same 
data but for W/t = 0.05. It is clear that as W/t is enhanced, the large density of states 
present at ~ Agp vanishes. 

III. FINITE DOPING 

We now consider finite doping. Here, we are confronted with a sign problem which leads, 
for a given accuracy, to an exponential increase of CPU time as the lattice size and inverse 
temperature grows. 

We first consider the particle number as a function of chemical potential at U/t = 4 
for W/t = 0,W/t = 0.05 and W/t = 0.15 on an 8 x 8 lattice in a temperature window 
T = 0.5t to T = 0.2t. The data is plotted in Fig. (^. Finite size effects grow as a function 
of W/t. This may be seen for example in Fig. where it becomes increasingly hard to 
extrapolate to the thermodynamic limit for large values of W/ 1. To give the reader a feeling 
of the size effects involved, we have plotted at W/t = 0.15 and T = 0.25t, some data for an 
L = 10 lattice (see Fig. ^ For the discussion of the data, we will omit the size effects in the 
considered temperature range. The dashed lines correspond to the value of the quasiparticle 
gap obtained in the zero-temperature and thermodynamic limits, (see Fig.jip. As discussed 
above, we expect for our model Agp = /ic. Hence, at T = 0, n{fi) = 1 for /i > /ic and 
n(/i) < 1 for fi < fic- From this construction, one may see that temperature effects are much 
more pronounced at W/t = than at W/t = 0.15. This result is easily understood within a 
single-particle picture mentioned in the introduction. The flatness of the dispersion relation 
signaled by the correlation length exponent u = 1/4 at W/t = leads to large temperature 
effects in n{fi). The transition to the description of the insulator-metal transition with 
exponent z/ = 1/2, leads to the suppression of large temperature effects in n(/i). At W/t = 
it was possible to pin down the exponent of the charge susceptibility, Xc ~ |/^ ~ f^cl'^^"^ only 
with the use of zero-temperature algorithms such as the PQMC . As argued above, at 
finite values of W, temperature effects are reduced. Nevertheless, it is still difficult to pin 
down the charge susceptibility exponent from finite temperature data. 
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We conclude this section by giving some numerical evidence, that the cause of the tran- 
sition in the nature of the insulator-metal transition is related to the enhancement of 
pairing correlations in the metallic state. The vertex contribution to the equal time pair-field 
correlations is given by: 

cr,S,S' 

(20) 

where 

P,,,(f) = {Al{r)A,M (21) 

with 

Here, fs0) = 1 and fd{S) = 1(— 1) for S = ±d^ (=tOy)- Per definition, P^si"^ = in the 
absence of interactions. Fig. |^ plots PJ.,(L/2, L/2) at (n) = 0.78, and U/t = 4 as a function 
of temperature. We consider an 8 x 8 lattice with W/t = 0, W/t = 0.05 and W/t = 0.15. 
As W/t grows, there is a substantial increase in PJ(L/2,L/2). This stands in contrast to 
the signal obtained in the extended s— wave channel where virtually no enhancement as a 
function of W/ 1 is seen. 



IV. DISCUSSION AND CONCLUSIONS 

We have given numerical evidence showing that the nature of the metal insulator tran- 
sition in the two-dimensional Hubbard model is radically altered by the inclusion of the 
W^— term. The low energy physics contained in this term was argued to be described by the 
three site terms obtained in a strong coupling expansion of the Hubbard model. The pre- 
sented data support the interpretation that the change in the nature of the doping induced 
metal insulator transition as a function of W/t is related to the onset of d— wave supercon- 
ductivity in the metallic state at finite W/t. The critical value of W, W^/t, at which the 
nature of the metal-insulator transition changes is uncertain. We can at present only give 
an upper-bound: W^ /t < 0.05 for U/t = 4. 

Our result are best described under the assumption of hyperscaling [II],^ where it is 
postulated that the singular part of the free energy scales as /^(A) = A'^^'^"'"^^and that there 
is a single relevant length scale ^. Here, A = |/i — fid is related to the length ^ ~ A"'", 
d is the dimensionality and z the dynamical exponent. Under the above assumptions, the 
charge susceptibility Xc, charge stiffness D, and localization length ^i, satisfy: 

e«~A-^ Xc-A^^'^-"), D ~ A''('^+"-2), (23) 

Since the control parameter. A, corresponds to the chemical potential, one obtains the 
additional constraint uz = 1 as well as (5 ~ ^ § being the doping concentration. 
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Thus, in the case of a doping induced metal-insulator transition, there is only one free 
parameter. The generic band metal-insulator transition does satisfy the above assumption 
and belongs to the z = 2, i/ = 1/2 universality class. We also note that the doping induced 
metal-superfiuid transition for two-dimensional bosons, is described by hyperscaling with 
exponents: v = 1/2 and z = 2 For the two-dimensional Hubbard model at f//t = 4, 
we have two numerical results, the localization length as well as the charge susceptibility, 
which are consistent with the above scaling relations provided that z/ = 1/4 and z = 4. At 
W/t = 0.05 and U/t = 4, our results for the localization length, ~ |yU — /^d"^^^ imply 
that u = 1/2. Furthermore, from our finite temperature results for it is plausible to 

assume that Xc ~ ^ f^c\^ the zero-temperature and thermodynamic limits as would be 
expected if the hyperscaling assumption is satisfied. For the confirmation the hyperscaling 
assumption, the calculation of the charge stiffness at T = as a function of doping and W/ 1 
is crucial. The prediction that D ^ at W = and that D 6 at W > remains for 
further studies. 

The above analysis suggests that the metal-insulator transition characterized by the 
dynamical exponent z = 4 is destabilized by a finite value of W to the superconductor- 
insulator transition with z = 2. When the dynamical exponent z changes from 2; = 4 to 
2; = 2 by inclusion of a finite amplitude of the W-term, the thus obtained wider dispersion 
relation causes a strong kinetic energy gain. This may be the origin of the instability of 
the incoherent metals with 2; = 4 to the superconducting state with z = 2. A flat band 
structure itself could be mimicked by a simple modification of the original band structure 
(e.g the inclusion of next-nearest neighbor transfer). However, to understand the instability 
due to the change in the universality class and the origin of the release of the suppressed 
coherence (i.e. D ~ (5^), it is crucial to notice that the z = 4 metals are a consequence of 
strong correlation effects and have to be constructed not from simple band electrons but 
from strongly renormalized charge excitations around the (±7r, 0) and (0, ivr) points in the 
Brillouin zone 

We have previously shown that at T = 0, (n) = 1, and U/t = 4:,W drives the system from 
an antiferromagnetic Mott insulator to a dx2-y2 superconductor. This quantum transition 
occurs at Wc/t ~ 0.3. In the superconducting state at half-band filling, charge and spin 
degrees of freedom have been studied and show rich physics [l^. The data shows the 
extreme compatibility between antiferromagnetic fluctuations and d-wave superconductivity. 
The relation between this state and the state obtained by doping the antiferromagnetic Mott 
insulator at U/t = 4, < W < Wc is intriguing. High-Tc cuprates are obtained by doping 
the Mott insulator. Some aspects of the experimental data seem to show the remnant 
of extreme compatibility between magnetism and c?— wave superconductivity. In nuclear 
magnetic resonance studies of YBa2Cu40s and HgBa2Ca2Cu30s+5, T2G of Cw-sites, which 
measures predominantly Rex(^ = (ti", tt), = 0) continues to grow monotonically with 
decreasing temperature even below the temperature scale at which 1/TiT starts decreasing 
25| . Although this feature is not clear cut in the single- layer-type mercury compounds |26 



the tendency of stronger growth of I/T2G than 1/TiT above seems to be universal. An 
extreme case of this feature is seen in our numerical calculations in the superconducting state 
at (n) = 1 where in the low temperature limit, 1/TiT scales to zero while Rex(^ = (tt, vr), ti; = 
0) diverges. The most remarkable feature is that the dynamical spin structure factor, S{q, u), 
at q = (tt, tt) has has a peak structure with diverging weight at finite frequencies . In 
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La2-xSrxCuOi^ S{q,uj), below Tc shows a sharp peak in momentum space centered at the 
incommensurate wave- vector, Qs = (vr, 7r)+5(7r, 0) and at relatively high energies: u > 7meV 
||27|| . This indicates that, the magnetic length scale at relatively high energies becomes large 
( > 50A ) as the temperature is lowered below T^. In YBa2Cu30-j_x, a sharp peak in 
frequency at 30 — 40 meV is observed in the superconducting phase p8| , p9| . These features 
again share some similarity with the t — U — W model at [n) = 1 in the superconducting 
state. At temperatures below the Kosterlitz-Thouless transition, spectral weight in S{q,uj) 
is transfered from low frequencies to higher frequencies. In the low temperature limit, 
S{q = (7r,7r),a;) shows a sharp peak in momentum space at finite frequencies. This peak is 
associated with the coherent part of the spin response since it is related to the long distance 
power-law decay of the equal time spin-spin correlations P^ . 

In conclusion, the t — U — W model at finite values of W/t appears to be a promising 
minimal model for high-Tc superconductors. Clearly, further studies for the complete un- 
derstanding of physical properties of the model as function of W, doping and temperature 
are required. 
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FIG. 1. The quasiparticle gap Agp = E^^^ — Eq as a function of inverse linear lattice length 
L at [/"/t = 4 and three values of W/t. Here, N = L^. The solid lines correspond to least square fits 
to the form a + h/L. The data points at l/L = 0, are the so extrapolated values of the quasiparticle 
gap. At W/t = 0.15 the lattice sizes L = 12, 10 and 8 were taken into account for the fit. 
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FIG. 2. The single particle Green function at half-filling R = (L/2, L/2), U/t = 4, W/t = 0.05 
as a function of imaginary time r. Due to particle-hole symmetry, and the fact that the considered 
values of L are even, G{R,—t) = —G{R,t). It is interesting to compare this data to the case 
W/t = 0(0]). At W/t = 0.05, an added electron propagates quicker in imaginary time than at 
W/t = 0. 
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U/t = 4, W/t = 0.05, (n) - 1 




FIG. 3. G{R = (L/2, L/2), a;) as obtained from equation |l^ as a function of L. Here we 
consider the parameter set U /t = 4, W/t = 0.05 values of w satisfying u) < Aqp ~ 0.62t The solid 
lines are least square fits to the form aexp f — L/-v/26)- 
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U/t = A 




FIG. 4. The localization length as a function of |// — ^d- The data at W/t = 0, U/t = 4 is 
reproduced from Ref. Q. The data at W/t = 0.05, C//t = 4 is obtained from Fig. (|3|). The solid 
lines correspond to the form |^ — ^d"*^ for u = 0.5 and v = 0.25. The data is consistent with 
l/i - Atcl"^/^ at VF/t = 0.05 and |/i - ^d"^/^ at W/t = 0. 
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N{u;), U/t = i, (n) = 1, T = 
W/t = 0.1b W/t = 0.0 




N{u;), Ult = A, (n) - 1, T-0 
W/t = 0.05 
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FIG. 5. T = density of states, N{uj) at half-band filling (i.e. = ) for = 4 and (a) to 
(e) W/t = 0.15, (f) to (j) W/t = and (k) to (n) W/t = 0.05. 
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FIG. 6. Number of particles per site, n, as a function of chemical potential for temperatures 
ranging from T = 0.2t to T = 0.5t. The dashed vertical line corresponds to the quasiparticle gap 
at r = and in the thermodynamic limit. This data is obtained from Fig. |l[ We consider an 
8x8 lattice, U /t = 4 and (a) W/t = 0, (b) W/t = 0.05 and (c) W/t = 0.15. In Fig. (c) we have 
included some data points for a 10 x 10 lattice. 



18 



U/t = 4, (n) = 0.78, L = 8 



Wlt = Qi 



V:/^(L/2,L/2) 
A:PJ(L/2,L/2) 



{a) 



Wit = 0.05 




(&) 



W/t = 0.15 




A ^ 



0.1 



FIG. 7. Vertex contribution to the pair-field correlations at the largest distance on an L = 8 
lattice in both the d— wave and extended s— wave channel. The data is plotted versus temperature. 
Here we fix the the particle number to (n) = 0.78, the Hubbard repulsion to U/t = A and consider 
(a) W/t = 0, (b) W/t = 0.05 and (c) W/t = 0.15 
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